Spontaneously axisymmetry breaking phase 
in a binary mixture of spinor Bose-Einstein condensates 
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We study the ground state phases for a mixture of two atomic spin-1 Bose-Einstein condensates 
(BECs) in the presence of a weak magnetic (B-) field. The ground state is found to contain a 
broken-axisymmetry (BA) phase due to competitions among intra- and inter-species spin exchange 
interactions and the linear Zeeman shifts. This is in contrast to the case of a single species spin- 
1 condensate, where the axisymmetry breaking results from competitions among the linear and 
quadratic Zeeman shifts and the intra-species ferromagnetic interaction. All other remaining ground 
state phases for the mixture are found to preserve axisymmetry. We further elaborate on the ground 
state phase diagram and calculate their Bogoliubov excitation spectra. For the BA phase, there exist 
three Goldstone modes attempting to restore the broken U(l) and SO(2) symmetries. 

PACS numbers: 03.75.Mn, 03.75.Kk, 67.60.Bc 



I. INTRODUCTION 



The spin degrees of freedom is described by the SU(2) 
symmetry group, which can take different irreducible rep- 
resentations, each corresponding to different values of the 
spin. For an atom in a hyperfine spin / state, its corre- 
sponding SU(2) representation is 2/+1 dimension. With 
the help of optical traps, spinor condensates with atoms 
in all these representations have been realized in many 
different atomic species and states 

Extensive experimental and theoretical investigations 
have focused mainly on atomic spinor condensates with 
spin equal to 1, 2, and 3. For the simplest case of 
spin-1 condensates 0-flll, the mean-field ground states 
% 12], and exact full quantum eigenstates [H, EH are 
intensively explored, also noted is the possibility of topo- 
logical fractionalized 1/2-vortex excitations in the anti- 
ferromagnetic phase [14| and the appearance and forma- 
tion of spin domains in the ferromagnetic phase (lH fl6l j . 
In addition, theoretical studies [§, [l7|, [3 either based 
on mean-field approximation or including quantum spin 
fluctuations, and experimental efforts [19i . |20| have been 
made to investigate the spin dynamics of the conden- 
sate. Similar efforts are directed at spin-2 condensates 
[HI, [2l], 22], whose ground states contain one extra phase: 
the cyclic phase, with an exotic symmetry described by 
the non-Abelian tetrahedral group T, potentially of sig- 
nificant value for demonstrating various topological ex- 
citations (23rl25j . Under the single spatial mode approx- 
imation (SMA), the exact eigenspectra and eigenstates 
of a spin-2 condensate are derived [12, [22J, which shows 
that the ground state is in general fragmented. The spin 
dynamics are also studied in the presence of a B-field, 
with the linear and quadratic Zeeman shifts [!, [2(| H?} ■ 
For the case of the spin-2 hyperfine manifold of a 87 Rb 
condensate 27], its spin dynamics show amplitude reso- 
nance when the B-field is tuned experimentally just like 
spin mixing in a spin-1 condensate [l?], EH- When the 



spin increases to 3 [6J, the ground state contains many 
more phases 0, [29j , and new possibilities for richer and 
more colorful physics become available, waiting to be fur- 
ther explored. 

Parallel to these studies of spinor condensates, the ef- 
forts are multi-component quantum gases have also re- 
ceived equal attention, most notably the case of pseudo 
spin- 1/2 condensates (30l435j j involving two atomic 
species or components. Since the total spin is conserved 
under the density dependent interaction between the two 
components, even in the presence of evaporative cooling 
if the interaction is SU(2) symmetric, A. B. Kuklov et al. 
[33| predicted that in the ground states all atoms con- 
dense into two orthogonal spatial orbitals, resulting in 
a maximally entangled many atom state. Replacing the 
two orbitals with two different species of atoms, Shi et 
al. [35| found a ground state condensate with an entan- 
gled order parameter. For two bosonic species of spin- 
1 atoms, such as 87 Rb and 23 Na atoms in the lowest 
hyperfine state manifold, the mixture reduces to a bi- 
nary mixture of spin-1 BECs. In the absence of external 
B-field, the ground state phases are essentially known, 
from the semi-classical phase diagram under the mean- 
field approximation [36| to the quantum results under the 
SMA for each of the two species and including quantum 
spin fluctuations 37-39]. Several interesting results have 
been revealed already. For example, the ground state 
reduces to a maximally entangled states [33, HI] where 
the mixture becomes fragmented and the number fluctu- 
ations [Hj] exhibit drastically different features from the 
isolated spin-1 condensates, when the inter-species' anti- 
ferromagnetic spin exchange interaction is large enough 
and when the inter-species singlet pairing interaction is 
ignored. 

In this work, we continue our study of a binary mixture 
of spin-1 condensates, in the presence of a weak external 
B-field. First, as in all experimental setups, a nonzero 
B-field inevitably exists no matter how small it may be. 
Second, it is motivated by the possibilities of interesting 
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physics due to the competitions between the linear and 
quadratic Zeeman shifts and the ferromagnetic interac- 
tions, which for the case of isolated spin-1 condensate 
is shown to induce a spontaneous axisymmetry breaking 
phase 0, 113] • We find that by tuning the inter-species 
spin exchange interaction and the linear Zeeman shift, 
and largely ignoring the quadratic Zeeman shifts (be- 
cause the B-field is weak), there exists an analogous BA 
phase in a system of binary mixture of spin-1 conden- 
sates. 

This article is organized as follows. In section II, we 
introduce the model formulation for a binary mixture 
of spin-1 condensates in the presence of an external B- 
field. Sec. Ill is devoted to the mean field results of the 
broken-axisymmetry (BA) phase, where we discuss the 
special case and the general case respectively. In Sec. 
IV, we elaborate on the Bogoliubov spectra for the vari- 
ous quantum phases of the mixture, and we point out the 
interesting zero energy Goldstone modes responsible for 
restoring the broken continuous symmetries: U(l) and 
SO(2). Also we propose a new scheme for classifying the 
various phases in the ground states using the number of 
nonzero components for their associated order parame- 
ters. Finally we conclude in Sec. V. The appendix pro- 
vides some mathematical details. 



II. OUR MODEL 



In this section, we first review the basic formulation in- 
troduced earlier for the inter-species atomic interactions. 
Between two distinguishable spin-1 atoms (of different 
species), they can be parameterized by the three diffcr- 

(12) 

ent scattering lengths a 1 2 which refer to three chan- 
nels of total spin F to t = 0, 1, and 2 respectively. Define 

9o% = 47rft2 aoli 2) 2 /A 1 with M = M 1 M 2 /(M 1 + M 2 ) the 
reduced mass and M; the mass for each atomic species 
(/ = 1,2), the corresponding pseudo-potential is given by 

F 12 ( ri - r 2 ) = (.g( 12 V + g[ 12) Vi + sf W(r x - r 2 )/2. 
^0,1,2 is the corresponding total spin projection operator, 
and Fx -F% = V2—V\ — 2Vq. The inter-species interaction 
is then expressed as 

Via(ri - r 2 ) = i(a + /3F X • F 2 + 1 V )5{r l - r 2 ), (1) 



where a = (g[ 12) + g£ 2) )/2, (3 = (_^ 12 ) + 5 <, 12 >)/2, and 
7 = (2 5 <, 12) - 3<?f 2) + <? 2 12) )/2- Denoting tf Mj ,(r) and 
$m f (r) as the two species annihilation field operator at a 
position r respectively, the Hamiltonian of a homogenous 



system of binary mixture under B-field becomes 
H = H1 + H2 + H12, 

Hi = Jdr 2^V 2 -pim + 9l m 2 )# m 

H 12 = - J dr {a#]$t$^ 

+&l& k F Uj • Faw^ + ^7 s f s| , (2) 

where pi and q\ are the linear and the quadratic Zee- 
man shifts, respectively, and summation over repeated 
indices are assumed. iT 2 has the same expression as Hi 
except for the substitution of subscript 1 by 2 and if! by 
4. s — (*i$_i - ^0*0 + is the inter-species 

singlet pairing operator. 



III. MEAN-FIELD GROUND STATES 

This study concerns the limit of weak external B-fields, 
when the quadratic Zeeman shifts can be neglected, or 
qi = g 2 = 0. Without loss of generality we assume 
P1P2 > 0. For the mixture of two alkali species, such 
as the spin-1 87 Rb and 23 Na condensates, under the as- 
sumption of interspecies interaction mainly arise from 
the contributions of the two valence electrons [4ll - l43| . 
we can approximate 7 = [44| . For a uniform sys- 
tem, it is convenient to expand the field operators in 
terms of plane waves as \& m = ^ k e ik r ak,m an d 

$ m = fi^ 1 / 2 ^ k e lk r 6k,m, where Q is the system vol- 
ume, and m (6k jTO ) denotes the annihilation operator 
of an atom of species one (two) with momentum Kk and 
magnetic quantum number m. Then the Hamiltonian of 
Eq. @ can be rewritten as 

k,m k 

+ 2iE :f bk f l,-k:, 

k 

k k 

where ei k = ft 2 k 2 /2Mi, pi )k = S q , m "q+k,m^q>™' and 

fi,k = Eq^^VicmftaV' Tne symbol : : repre- 
sents the normal ordering of the operators. Similar to 
Eq. @, H2 is given by replacing the subscript 1 by 2, 
and the operator a by b in Hi . 

The condensate component corresponds to the zero 
momentum k = state, which is occupied by a macro- 
scopic number of atoms. The corresponding Hamiltonian 
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becomes 



-ffBEC = / ] -m{jp\OQ , m ao,m + P2feo, m &0,m) 
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Under the mean-field approximation, the operator oo, m 

(Vm) is replaced by c-number \/A?iCi,m {VN 2 C,-z,m), 
where N\, iV 2 are the atom numbers for the species one 
and two respectively, and Ci = (Ci,1j Ci,0> Ci,-i) T and 
C2 = (C2.i1 C2.0, C2,-i) T are normalized spin-1 spinors. 
As a result, in the mean-field approximation the ground 
state spinor wavefunction is found by minimizing the 
mean-field spin-dependent energy given by 



En = 



+ ^2n 2 ^(F 2 ) 2 + i/3n 2 (F 1 ).(F 2 ) 

= -pfu - xpf 2z + I (f3[ f 2 + p' 2 ff + p% ■ h), (5) 

where nj — Nj/Cl is the condensate density of the j- 
th species and fi, 2 = (Fi, 2 ), p\ = V2 N 2 /{N 1 x) = p, 
p[ = /Sim, f3' 2 = famN-i/Nx and p' = /3n 2 . The pro- 
cedure to find the mean-field ground states is summa- 
rized in the appendix^ where we discuss the structure 
of the ground states classified into two classes: with or 
without broken-axisymmetry. According to Ref. a 
state with a nonzero transverse magnetization is called 
a broken-axisymmetry phase. Analogously states pre- 
serving axisymmetry refer to those with zero transverse 
magnetization. In the BA phase, both species are found 
to be fully polarized. For axisymmetry preserved ground 
states, the wave functions take the form 



e^V( 1 + /j Z )/ 2 


V(Wjz)/2 



(6) 



where Xj and <fj are arbitrary phase angles, and j = 1, 2. 
When /i z = or / 2z = 0, the ground states include an 
infinite family of degenerate states with that in Eq. (|6|). 
More rigorously, we note that the nematic order for the 
above state (J6j) actually is not axisymmetric For the 
BA phase, the ground state wave functions take the form 



/ 



e -i Vj cog 2 & 

•\/2 cos sin 
e 1 ^ sin 2 ^ 



(7) 



where Xj and fj are arbitrary phase angles, and (p 2 — 
ip 1 +7r(mod 2?r), j = 1,2. 9\ and 2 are determined by f± z 
and / 2z respectively as cos Q\ — j\% and cos 8 2 = f 2z . The 
ground state phase diagram is characterized into three 
cases according to the intra-species spin exchange inter- 
action parameter P[ and P' 2 . Without loss of generality, 
in the following we assume \P[ \ < \P' 2 \. 
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FIG. 1: (Color online). The ground state phase diagram of 
our model system of a binary spin-1 condensate mixture at 
fixed values of f3[ and P' 2 for the special case of x = 1. The 
black dash-dotted lines correspond to p = 0, which serve as 
guides for the eye. Blue (Red) lines denote continuous (dis- 
continuous) phase boundaries. The first three subplots denote 
fixed intra-specie spin exchange interaction parameters at 
&)/Wi\ =■ (a) (-1,-2); (b) (1,2); and (c) (-1,2). The last 
one (d) illustrates the dependence of ground state order pa- 
rameters /i z (blue solid line) and /2 Z (red dashed line) of the 
BA phase on the linear Zeeman shift p at a fixed value of 
P' = 7\P[\. 



A. the special case of x = 1 

In this subsection, we consider the special case of x = 
1, where the parameter defined as the partical number 
multiplying the linear Zeeman shift, of the two species are 
equal (N1P1 = N 2 p 2 ). The ground states are as shown in 
the Fig. ([1]), where the first three subplots denote fixed 
intra-specie spin exchange interaction parameters at (P[ , 
P'2)/\P'i\ =■ (a) (-T-2); (b) (1,2); and (c) (-1,2). The 
last one (d) illustrates the dependence of ground state 
order parameters /i z and / 2z of the BA phase on the 
linear Zeeman shift p at a fixed value of P' — 7\P[\. 

First, we discuss the case when the two spin-1 conden- 
sates are both ferromagnetic (P[ < and P' 2 < 0). The 
ground state has two phases: the FF phase and the BA 
phase. The FF phase is a typical axisymmetry preserv- 
ing phase with /i z = /2z = sign(p), which is the same 
as that found in [36| with atomic spins in each species 
fully polarized and aligned parallel to each other with 
ff = ff — fi ■ f 2 = 1. In the BA phase, the two spin vec- 
tors of each species are remain fully polarized but now 
tilted at an angle 0\ and 2 with 61=62 = arccos(p//3') 
with respect to the z-axis. This is illustrated in the Fig. [5] 
From Eq. (|A3|) . we confirm that in this case the BA phase 
exists only when /3' > 0, which is consistent with numer- 
ical and analytical results shown in the Fig. [Ha). The 
boundary between the FF phase and the BA phase can 
be derived by comparing the energy of Eq. ([5]), which 



FIG. 2: (Color online). A schematic illustration of the 
spin vectors for both condensate species in the broken- 
axisymmetry phase. The weak external B-field is along the 
z-axis direction. The blue (red) arrow denotes the spin vec- 
tor of species one (two), which tilts at an angle with respect 
to the z-axis. The two spin vectors and the z-axis are in the 
same plane. The blue and red dotted dash lines correspond to 
the orbits of the end points of two spin vectors in the degener- 
ate ground state of the broken-axisymmetry phase. The black 
dotted dash lines represent the transverse magnetizations f 1 ± 
and f2± of two species respectively. 

turns out to be \p\ — (3' . In Fig. Hfd), we present the 
dependence of the order parameters f± z and / 2z on the 
linear Zeeman shift p at a fixed value of j3' = 7|/3jJ. In 
the region of \p\ < /?', /i z and / 2z change smoothly from 
— 1 to 1 simultaneously. 

We then turn to another case when one of the spin-1 
condensates is ferromagnetic (f3[ < 0) and the other is 
antifcrromagnetic or polar {j3' 2 > 0). In the absence of B- 
field, the competition between the intra- and inter-species 
spin exchange interactions gives rise to three phases [36[ . 
When /3 < — 2f3 2 , the ground state is in the FF phase. 
Increasing the value of /?, there arises another phase, 
where the spin vector in the ferromagnetic one is fully 
polarized, and the other in the polar one is not fully 
polarized, giving rise to ff = 1, f| = ((3' /2(3 2 ) 2 , and 
f i • f 2 = -/3/2/3 2 . In the other limit when /?' > 2/3 2 , 
the ground state is in the AA phase where the two spin 
vectors are fully polarized into opposite directions with 
fp = ff = 1 and fi • f 2 = — 1. In the presence of a weak 
external B-field, these three phases for the ground state 
in the phase space will extend into three regions. The 
first one is the FF phase, which has the same properties 
as that in the p — case, except now all spins are polar- 
ized along the external B-field direction. The second one 
corresponds to the other axisymmetry preserved phase 
denoted by P+N, where we use symbol P and N to rep- 
resent respectively the fully and non-fully polarized spin 
vectors of two species. The subscript + is used to clarify 
that the fully polarized spin is along the same direction 
of the B-field. The other remaining one is the BA phase, 
which is totally the same as that in the case of two ferro- 
magnetic spin-1 condensates. In Fig. ffljc), we illustrate 
the ground state phase diagram. The boundary between 
the FF phase and the BA phase is the same as that in 




P'flP'J) PflP',1) 
FIG. 3: (Color online). The same as in Fig.[T] but for x = 2/3. 



Fig. QJa), except that when 0' < 2(3 2l the ground state 
is changed to the P + N phase, which also means that the 
boundary between the P + N phase and the BA phase is 
P' = 2/3 2 , which is consistent with Eq. (|A3|) . Only when 
P' > max(0, 2/3' 1 x, 2/3' 2 /x) = 2j3' 2 is satisfied will the BA 
phase possibly exist. The remaining boundary between 
the FF phase and the P + N phase satisfies the relationship 
\p\ — (3' 2 + /3' /2. In the P + N phase, we have /i z = sign(p) 
and/ 2z = (p-/ lz /?'/2)//3 2 . 

Next, we consider the last case with two spin-1 po- 
lar condensates ((3' 2 > (3[ > 0). In Fig. QJb), we 
present the phase diagram as well as the boundaries 
separating the phases. Similar to the second case of 
Fig- Die), the ground state has three phases: FF, 
P+N, and BA. The order parameters and the bound- 
aries between these three phases are the same as that 
in the previous case of Fig. [ljc). The only difference 
arises from the competition with the intra-species in- 
teraction in the anti-ferromagnetic condensate. A new 
phase emerges in the center of the P + N phase, which 
is denoted as the NN phase. In this new phase, nei- 
ther spin vectors of each species are fully polarized. 
The boundaries between the NN phase and the P + N 
phase satisfy the following relationships ((3' /2 — (3 2 )\p\ — 
l3' 2 /4 - p[p 2 . For the ground states in the NN phase, 
we have / lz = (0'/2 - &)p/ \f3' 2 /4 - and / 2z = 

{p/2-K)p/(j¥ 3 H-F 1 &). 



B. the general case 

The ground state phase diagram for the special case 
of x = 1 as discussed above contains both axisymmetry 
preserved and broken phases. Their boundaries can be 
partially determined by Eq. (|A3|) . The AA phase is found 
to exist only when p = 0. It amounts to an infinite set of 
degenerate states, where the two spin vectors are antipar- 
allel while remaining fully polarized. As long as x ^ 1 is 
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FIG. 4: (Color online). The same as in Fig.Q] but for x = 3/2. 



satisfied, the region of the AA phase extends to p ^ 0, 
where one spin vector will be parallel to the B-field, while 
the other one become antiparallel. The order parame- 
ters (/iz,/2z) for the P+N, NN and BA phases change 
into the following forms (sign(p), (2xp~ sign(p)fl')/2(3' 2 ), 
((2xpp> - 4p(3' 2 )/W' 2 - W 2 ), (2pf3> - 4xp(3[)/((3' 2 - 
4/3^)) , (xp//3>- (x 2 - 1)0' /4xp, p/(3' + (x 2 - l)(3'/4x 2 p) , 
respectively. 

We first consider the case with x < 1, whose ground 
state phase diagram is demonstrated in Fig. [31 Due to 
our choice of \/3[\ < \j3' 2 \, the ground state phase diagram 
is similar to the case of x = 1. One notable difference 
is that the AA phase appears in the middle of the BA 
phase, and the boundaries between them are given by 
\p\ = — x)/2x. And the boundary between the BA 
phase and FF phase is also changed to be described by 
\p\ = j3'(x + l)/2x. As a result of the appearance of 
the AA phase, the P + N phase and the A A phase will 
now in touch with each other, and their boundary is 
given by |p| = (/?' — 2f3' 2 )/2x. The NN phase remains 
immersed in the P + N phase, and their boundary sat- 
isfies |p| = (/3' 2 - 4:/3[l3 2 )/(2x(3' - 4/%). In Fig. [3J we 
illustrate the phase diagram, where the blue (red) solid 
lines denote the continuous (discontinuous) phase bound- 
aries. In Fig.[3Jd), we display the order parameters of the 
ground states at fixed interspecies spin-exchange interac- 
tion j3' = 7\fi[\. Due to the emergence of the AA phase in 
the middle of the BA phase, f 2z changes from -1 to 1 twice 
when p changes from — oo to oo, while /i z oscillates near 
the maximum and minimum values respectively, which is 
different from the special case of x discussed previously, 
where /i z and f 2z change simultaneously from -1 to 1. 

Next we discuss the other case of x > 1 , whose ground 
state phase diagram is demonstrated in Fig. SJ Due to 
the choice of \f3[\ < \f3' 2 \ and x > 1, the ground state 



phase diagram is much more complicated than that of 
x < 1. For two ferromagnetic condensates (f3[ < and 
P 2 < 0), the ground state phase diagram is similar to 
that of x < 1. With one ferromagnetic (/3J < 0) and one 
anti- ferromagnetic condensate (f3' 2 > 0), the phase dia- 
gram becomes different from the case of x < 1. A new 
phase denoted by P_N emerges, where the subscript — 
is used to denote a spin vector being fully anti-aligned 
with the external B-field. This phase arises from in- 
creased interspecies spin-exchange interaction j3' in this 
case of x > 1. Starting from the P + N phase, with the 
increase of /3' , the spin vector of the second species be- 
comes incresingly polarized antiparallel to the B-field, 
which weights in more linear Zeeman energy contribution 
from the second species, and the spins for both species 
will flip to the opposite directions. The order parameters 
for the P_N phase take the forms /i z = — sign(p) and 
/2z = (2^p + sign(p)/3')/2/?2. Incidently, we find that the 
phase boundaries between the P + N phase and the BA or 
the P_N phases are the same: described by /?' = 2(3 2 /x. 
The remaining phase boundaries take the same forms as 
that in the case of x < 1, except for the one between 
the phase of P_N and AA, which takes the form \p\ = 
(2f3' 2 — f3')/2x. For the choice of two anti- ferromagnetic 
condensates (f3[ > and f3' 2 > 0), in addition to the P_N 
phase, a new phase denoted by NP + emerges between the 
BA phase and the NN phase. The order parameters for 
the NP + phase are given by /i z = (2p — sign(p)(3')/2(3' 1 
and f 2z = sign(p). The boundaries between the NP + 
phase and the NN, the AA, and the FF phases take the 
forma \p\ = (/3' 2 - 4^/3^/(2/3' - 4a^), \p\ = /3'/2 - ft, 
and \p\ — (3' /2 + f3[, respectively. The complete phase 
diagram is illustrated in the Fig.Q] Additionally we illus- 
trate the order parameters for the ground state at fixed 
value of /3' = 7|/3JJ. When x > 1, /i z is observed to 
change from -1 to 1 twice as p changes from — oo to oo, 
while f 2z oscillates near the maximum and minimum val- 
ues, in contrast to the special case of x < 1 as shown in 
the Fig.^d). 



IV. BOGOLIUBOV SPECTRA AND 
GOLDSTONE MODES 



This section is denoted to the discussion of Bogoli- 
ubov excitations of the ground state phases we derived 
above. We follow the numb er- conserving approach ini- 
tially developed by Gardiner et. al. [45( , where clq m 



A,™) 1/2 C 



(bo, m ) is replaced by (Ni - 

((N2 - J2 k ^ m bi m K m ) 1/2 (2, m ) |4Q|,|46|,|47|, and keep- 
ing only quadratic terms of the operators to derive the ef- 
fective Hamiltonian by using the same notations as Refs. 

Em, 
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Hf = (eik-pim+pi^iz) -/3i?ii(Fi) 2 )a ] j < ., m ak,m + A«i(Fi) • ^ 4,m F i™«a k ,„ 

k^O,m k^O, m,n 

+ Y^^l.yP^ + ^lA-k + ^I )k ^l,-k) + ^y 1 E( 2 ^lk • ^i,k + A,k ■ A,-k + ^.k • ^i,-k), 

k^O k^O 

#if = f E -( F i)-( F 2)(«2a^ m a k , m + 7 11 6j iim 6 k , m ) + ^i(F 1 ). ]T &k,m F 2mA,n 

k^0,m k^O, m,n 

+ ^( F 2) • E 4 )TO F lmnfik,n + E^U^2,k + ©l,k% + ©l,k©2,-k + ^I, k ^L k ) 

k/O, m,n k^O 

+ £(^?,k ' ^2,k + Ji,k • + A,k • + ^k • ^,_ k ) + £b, (8) 

k#0 



where £q is a constant term and we have defined £>i ;k = 

m n Cl,m F lm«^k,n' As before 
H 2 S takes the same form as Hf s except for the replacing 
of the subscript 1 by 2 and the operator a by b. 

Following the notation and approach of Ref. [40|, we 
define operators 

P k = («k4, <2k,0) Ok,-li &k,lj £>k,0) &k,-l) T , 

p k - (4,i.4 > o»4,-i ) &k,i^k > o.^,-i) r , (9) 

which satisfy the Heisenberg equation 

ih^-P k =M(k)K+M(k)P*_ k , (10) 

to calculate low-lying Bogoliubov modes. A4(k) and 
Af(k) are chosen as real with suitable real valued Ci and 
Q 2 and symmetric 6x6 matrix. The operators for Bo- 
goliubov quasi-particles are defined accordingly 

Qk = (Qk,l; 2k,2, Qk,3, Qk,4, Qk,5j Qk : 6) T 

= W(fc)P k + V(k)P*_ k , (11) 
which enforces a diagonalized effective Hamiltonian 

H cS = E EmQl^Qk^ + E vac , (12) 

k ,m 

where E v ^ c is the vacuum state energy for the quasi- 
particles. 

As the quasi-particle operators Q k>m and satisfy Bose 
commutation relations, we obtain 

«^Qk = £(k)Q k , (13) 

where E(k) is a diagonal 6x6 matrix: diag(-E) = 
(Ei , E 2 , E 3 , Ei, E 5 , E 6 ) T . The eigenvalues of the quasi- 
particles are found from solving the following eigenequa- 
tion [Hi 

(M+N)(M-N)(U + Vf = (U + V) T E 2 . (14) 



I 

We choose to reexamine the ground states of the 
Hamiltonian of Eq. (0| before discussing the Bogoli- 
ubov excitations and the Goldstone modes. As shown 
in the appendix [A] there exist two types of mean field 
phases in the ground state, one preserves axisymmetry, 
and the other breaks axisymmetry. The axisymmetry 
preserved phases include FF, AA, P+N, P_N, NP+, 
and NN phases, etc. All ground state phases can be 
reclassified into four classes according to the numbers 
of nonzero component in their order parameters: the 
FF and AA phases falls into the two-component class, 
the P+N, P_N, and NP + phases belongs to the three- 
component class, while the four-component class includes 
the NN phase and the six component class includes the 
BA phase. 



According to the Goldstone theorem: a gapless mode 
emerges whenever a continuous symmetry is sponta- 
neously broken. For the two-component class, two U(l) 
symmetries for the two atomic species are broken |49j. 
One and two more spin rotation symmetries are bro- 
ken in both the three-component and four-component 
classes. For the BA phase, two U(l) symmetries and 
one SO (2) symmetry are broken [501 ] - To verify the 
above arguments, we solve the eigen-equations of Eq. 
(1141) employing both analytic or numerical techniques. 
In the ground state of the Hamiltonian of Eq. (Q}, 
unequal densities of two species can modify the spin- 
exchange interaction and linear Zeeman shifts. With- 
out loss of generality, in the following numerical calcu- 
lation, we will assume Ni — N 2 and ni = n 2 = n. Be- 
cause spin-independent interactions do not change the 
general properties of the Bogoliubov excitations, we will 
set intra- and inter-species spin-independent interaction 
as OLiUi = a 2 n 2 — a^/nin 2 = an — 20\f3[ | and we further 
assume M% = 3Mi. 
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A. two-component class consider the FF phase by assuming Ci — C2 — (1, 0, 0) J 

and p > 0. After detailed calculation, we find 

As summarized above, the two-component class con- 
tains two phases: the FF phase and A A phase. We first 



El = ~ (e 2 k + e* k -yf (e 2 k - 4J 2 + 4(on + fin) 2 e lk e 2k ) , 

E l = \ (e? k + £ 2k + ]f (e? k - £ L) 2 + 4(«n + /3n) 2 £ lk £ 2k ) , 

1 / \ 2 

£3 = - (em + e 2 k - fin + pi + p 2 - \/ (eik - e 2k + pi - P2) 2 + (/3«) 2 J , 

2 / \ 2 

El = - (eik + e 2k - fin + pi + p 2 + \/ (eik - e 2k + Pi - P2) 2 + (l3n) 2 j , 

El = (eik - 2/3m - /3n + 2 Pl ) 2 , 

£ 6 2 = (e 2k - 2/3 2 n - /3n + 2p 2 ) 2 , (15) 



where e 2 k = e 2 k + 2(am + /3in)ei k and ef, k = e\\ + 
1{a 2 n + /3 2 n)e 2k . From Eq. ([T5"j) . we confirm that there 
are two gapless Goldstone modes with eigenvalues £1 
and £ 2 , associated with the two coupled U(l) symmetry 
breaking modes between the two Mp = 1 components of 
two species respectively. The two coupled modes between 
the two Mp = components are gapped, and their exci- 
tation spectra are denoted by £3 and £4. The remaining 
two excitation spectra are the two modes coupling be- 



tween two Mp = — 1 components, which are gapped as 
well. In the boundary between the FF phase and the 
BA phase, where \p\ — fin, we find the excitation spectra 
of £3 reduce to a free particle form at small values of 
momentum k. 

For the AA phase, without loss of generality, we can 
assume the order parameters of the forms £1 = (1, 0, 0) T 
and (2 = (0,0, 1) T . The excitation spectra can be ana- 
lytically retrieved in this case as 



E l = \ ( £ ik + 4< - \j (e? k - e L) 2 + 4 ("« - /3") 2 £ik£2k) , 

E l = g ( e ?k + 4< + \] (eik - e 2k) 2 + 4(an - fin) 2 e lk e 2k ) , 

El = - (eik - £2k + Pi + P2 - V (eik + &2k + fin + Pi - P2) 2 - {fin) 2 

El = - (eik - e 2 k + P\ + P2 + V ( £ ik + ^2k + fin + pi - p 2 ) 2 - (fin) 2 

El = (e lk ~2fi 1 n + fin + 2p 1 ) 2 , 

El = (£ 2k -2/3 2 n + /3n-2p 2 ) 2 . 



Again there are two gapless Goldstone modes, associated 
with the breaking of the two coupled U(l) symmetries 
between the components of (^1 and C2,-i denoted by E\ 
and £ 2 . The two coupled gapped modes of the Mp = 
spin components for the two species are denoted as £3 
and £4. The remaining two gapped modes come from 
the coupling between the components of £1.-1 and £2,1- 
All results are consistent with our previous discussions. 



(16) 

I 

B. three-component class 

For the three-component class, the spin vector of one 
species is fully polarized while the other one is partially 
polarized, resulting in the breaking of two U(l) symme- 
tries and one SO (2) spin rotation symmetry. For all 
phases: P+N, P_N, and NP + , their Bogoliubov spec- 
tra contain three gapless Goldstone modes. We verify 
the above results numerically by diagonalizing Eq. (|14l) . 
The case of the P + N is shown for illustrative purposes 
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FIG. 5: (Color online). The Bogoliubov excitation spectra for 
the P+N phase. The intra- and inter-species spin-dependent 
interaction parameters used are /3[ < 0, /3' 2 = 2|/3l|, and ft' — 
3\fti\. The B-field parameters are pi = p 2 = 2.5\ft[\. The 
inset shows the small momentum k region. 



0.45 



0.30 



CM E 

LU 



0.15 



0.00 




0.00 0.01 0.02 0.03 0.04 0.05 



25 
20 
15 



0.0 



0.1 



0.2 



0.3 



0.4 



FIG. 7: (Color online). The same as in Fig. [5] but for the 
BA phase. The parameters adopted are ft[ > 0, ft' 2 = 2\ft'i\, 
ft' = 6\fti\, and pi =p 2 = 4ft'/5. 



D. six-component class 

Finally for the six component class, the axisymmetry is 
broken and its order parameter taking the form of Eq. ([TJ . 
Due to breaking of two U(l) and one SO(2) symmetries, 
the Bogoliubov excitation spectra contain three gapless 
Goldstone modes, as confirmed by numerically solving 
Eq. (|14l) . xi — Xi — fi — and f 2 = ^ are used for 
numerical calculations, together with the intra- and inter- 
species spin-dependent parameters ft[ > 0, ft' 2 = 2 1 | , 
and ft' — Q\ft[\- The linear Zeeman shifts are assumed 
to satisfy p\ = p 2 = 4/3'/5. In Fig. [71 we present all 
Bogoliubov excitation spectra, from which one can easily 
verify the three gapless Goldstone modes. 



FIG. 6: (Color online). The same as in Fig. [5] but for the 
NN phase. The parameters adopted are ft[ > 0, ft 2 = 2|/3j|, 
ft' = \ft'x\, and px= P 2 = 1/3x1/2. 



in Fig. [5j where we assume spin-dependent parameters 
ft[ < 0, ft' 2 = 2\ft[\, ft' = 3\ft[\, and linear Zeeman shifts 
Pi = P2 = 2.5\ft[\. The order parameters take the forms 
of Eq. (j6]), where now we have /i z = 1 and we choose 

Xl = X2 = fl = f>2 = 0. 



C. four-component class 

There exists only a single one phase: the NN phase, 
for the four-component class. Its order parameter takes 
the form of Eq. ([6]). For simplicity, we assume xi = 
Xi = fi — f>2 — 0. In Fig. [6] we illustrates its collec- 
tive excitation spectra, where the intra- and inter-species 
spin-dependent interactions are parameterized by ft' > 0, 
02 = 2\ft[\, and ft' = \0[\. The B-field parameters are 
chosen as pi = p 2 — \0[\/2. Indeed one can infer that 
there are four gapless Goldstone modes in Fig. [51 consis- 
tent with the previous discussions. 



V. CONCLUSIONS 

In conclusion, we construct the ground state phase di- 
agram for a binary mixture of two spin-1 condensates 
under an external weak B-field. Our results show that 
the phase diagram can be classified into two categories: 
the axisymmetry preserved and broken phases. The first 
one is shown to contain the FF, AA, P+N, P_N, NP + , 
and NN phases, while the second one contains only sin- 
gle phase: the BA phase, which appears only when 
the inter-species anti-ferromagnetic spin-exchange inter- 
action is sufficiently large when compared to the intra- 
species spin-exchange and linear Zeeman shifts. Further- 
more, according to the number of the gapless Goldstone 
modes, the ground state phases can be classified into four 
classes: two-, three-, four-, and six-component class, re- 
spectively accompanied by two, three, four, and three 
gapless Goldstone modes. We calculate all the Bogoli- 
ubov excitation spectra with both analytical or numerical 
techniques, and confirm our arguments given above. 

Since the interspecies s-wave interaction parameters 
between two spin-1 atoms, such as between 87 Rb and 
23 Na atoms, are not precisely known, we cannot con- 
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elude which of the above discussed phases can be realized 
experimentally within presently available experimental 
technologies. Given the promising innovation of tuning 
interspecies interaction parameters with optical Feshbach 
resonances and other tools, it seems that phases we study 
can stimulate sufficient experimental interests to confirm 
their existence. 
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Appendix A: constructing the ground state phases 

In this appendix, we give the detail discussion on how 
to clarify all possible ground state phases for a binary 
mixture of spin-1 condensates under an external weak 
B-field. 

First, we assume that spontaneous axisymmetry break- 
ing doest not occur in either species, which results in 
fi± = f%± = 0. For species one, we then have 



0, 



(Al) 



+ d2 



0. For the latter 



Eq. (HH): Co = or Ci 

case, we infer that |Ci 1 ' ) | 2 — IC-il 2 an d /iz = 0, which is 
a state degenerate with Ci = (e iipi ,0,e iv ~ 1 ) T /V2, with 
arbitrary angles (pi and ip-i. Therefore the ground states 

is limited to the subspace with Cg 1 " 1 = 0. For the sec- 
ond species, an analogous reduction gives a similar sub- 

(2) 

space with Q =0. Summarizing, without spontaneous 
axisymmetry breaking, the ground states of the energy 
functional for a binary mixture of two spin-1 conden- 
sates, Eq. ([5]), can be chosen as Ci 

c 2 = (Ci 2 



(CP,0,C-l) T and 



(1)\T 



Next, we consider the structure of all possible sponta- 
neous axisymmetry breaking ground states. For a single 
spin-1 condensate, the linear Zeeman shift does not in- 
duce BA phase. If a spontaneous BA phase arises in 
a mixture, both species must break axisymmetry. We 
therefore search for the ground states of the BA phase 
in the subspace of ff = a 2 and f| = b 2 , with a 2 < 1 
and b 2 < 1. Since pip 2 > 0, the BA phase appears for 
anti-ferromagnetic interspecies spin-exchange interaction 
with P' > 0. In the ground states, the transverse mag- 
netization of the two species must align along opposite 
directions. 

The order parameter for the BA phase can be found 
from the following equations 



d£. 



d£ s 

df2z 



f 2z = 0. (A2) 



The solution of the above equations for the BA phase is 
given by /i z = xp//3' + a 2 ft' /Axp — xb 2 j3' /Ap and / 2z = 
p/(3' + b 2 (3'/Ap - a 2 f3'/Ax 2 p. We can calculate the hrst 
derivatives of the spin-dependent energy function £ s with 
respect to a 2 and b 2 at the BA phase, and the results are 



d£ s 
da 2 



for the ground state. Without loss of generality, an uni- 
tary transformation can be affected to set (q 1 ' as a real 
number. As a result, there left only two solutions for 



-A 



2z 



2^2^71 



1 



(2(3[ - p'/x), 



d£ s 
dP 



2^ 2 2^2^ 



flz 



f 2 

J 2z 



P'x). (A3) 



We therefore conclude that the BA phase is possiblely 
exist only when (3' > max(0, 2(3[x, 2f3' 2 /x). Because the 
first derivatives are both independent of a 2 and b 2 , the 
ground states of the BA phase will have a 2 — b 2 — 1. 

The actual ground state is found through a compari- 
son of the ground state energies for two types of phases: 
with/without broken- axisymmetry In addition, the re- 
sults found here are subsequently affirmed with the nu- 
merical ones from the use of simulated annealing method 
[36j to minimize the spin-dependent energy functional. 
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